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Radiation and MHD effects on an oscillatory flow of 

a Jeffrey fluid in a circular tube 
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Abstract — In this paper, we studied the effects of heat transfer 
and MHD on oscillatory flow of Jeffrey fluid in a circular tube 
with radiation. The expressions for the velocity field and 
temperature field are obtained analytically. The effects of 
various emerging parameters on the velocity field and 
temperature field studied in detail with the help of graphs. 

Index Terms — Jeffrey fluid, Reynolds number, Prandtl number, 
Grashof number, Nusselt number. 

I. Introduction 

The effect of heat transfer on unsteady MHD oscillatory flow 
of fluid in a tube is encountered in a wide range of engineering 
and industrial applications such as molten iron flow, recovery 
extraction of crude oil, geothermal systems. Many chemical 
engineering processes like metallurgical and polymer 
extrusion processes involve cooling of a molten liquid being 
stretched in a cooling system. Some polymers fluids like 
polyethylene oxide and polysobutylene solutions in a cetane, 
having better electromagnetic properties are normally used as 
cooling liquid as their flow can be regulated by external 
magnetic fields in order to improve the quality of the final 
product. Also, the radiative heat transfers is an important 
factor of thermodynamics of very high temperature systems 
such as electric furnaces, solar collectors, storage of nuclear 
wastes packed bed catalytic reactors, satellites, steel rolling, 
cryogenic engineering etc. The study of such flow under the 
influence of magnetic field and heat transfer has attracted the 
interest of many investigators and researchers. Unsteady and 
oscillatory flow of viscous fluids in locally constricted, rigid, 
axisymmetric tubes at low Reynolds number has been studied 
by Ramachandra Rao and Devanathan (1973), Hall(1974), 
Chaturani and Upadyaya (1978) and Schneck and Ostrach 
(1975). Haidar (1987) have considered the oscillatory flow of 
a blood through an artery with a mild constriction. Ogulu and 
Abbey (2005) have studied the effect of heat transfer on the 
motion of blood in a diseased artery has been modeled under 
the optically thin fluid assumption. Several other workers, 
Misra and Singh (1987), Ogulu and Alabraba (1992), Tay and 
Ogulu (1997) and Elshahed (2003), to mention but a few, 
have in one way or the other modeled and studied the flow of 
blood through a rigid tube under the influence of pulsatile 
pressure gradient. Arunakumari et al. (2012) investigated the 
influence of heat transfer on MHD oscillatory flow of Jeffery 
fluid in a Channel. Unsteady Flow of a Jeffrey Fluid in an 
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elastic tube with stenosis was investigated by Sreenadh et al 

( 2012 ). 

In view of these, we studied the effects of heat 
transfer and MHD on oscillatory flow of Jeffrey fluid in a 
circular tube with radiation. The expressions for the velocity 
field and temperature field are obtained analytically. The 
effects of various emerging parameters on the velocity field 
and temperature field studied in detail with the help of graphs. 


II. .MATHEMATICAF FORMUFATION 

We consider an oscillatory flow of a Jeffrey fluid in a heated 
uniform cylindrical tube of constant radius R in the presence 
of uniform magnetic field B 0 . The wall of the tube is 

maintained at a temperature T . We choose the cylindrical 

coordinates (V,(9, z) such that r = 0 is the axis of 

symmetry. The flow is considered as axially symmetric and 
fully developed. The geometry of the flow is shown in Fig. 1. 



mum 

0 

Fig. 1 The physical model of the 


Fig. 1 The physical model of the problem 
The equations governing the flow (Kavitha and 
Ramakrishna Prasad, 2012) are given by 


dw dp ju 


f dZ 1 


m = 


—— + 


dz 1 + 4 


8 w l dw 

+- 


A 


v 


dr 


J 


P c > 


dT_ 

dt 


= k 


o 


f 


V 


r dr 

( 1 ) 

d 2 T 1 dT 

- 2~ ^- 

dr r dr 


-aBl W +p g p{T-T f ) 


\ 


J 


dq 

dr 


( 2 ) 


where p is the fluid density, is the ratio of relaxation to 
retardation times, /J is the fluid viscosity, p is the pressure, 
W is the velocity component in Z - direction, g is the 

acceleration due to gravity, C is electrical oconductivity, j3 
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coefficient of thermal expansion, T is the temperature, k 0 is 


the thermal conductivity and C p is the specific heat at 

constant pressure. Following Cogley et al. (1968), it is 
assumed that the fluid is optically thin with a relatively low 
density and the radiative heat flux is given by 

dq 


dy 


= 4a 1 (T, J ,-T) 


(3) 


here CC is the mean radiation absorption coefficient. 

The appropriate boundary conditions are 

w — 0, T — T w at r = R 
dw 

— = 0, T = T x at r = 0 (4) 

or 

Introducing the following non-dimensional variables 
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into the Eqs. (1) and (2), we get (after dropping bars) 
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where Pr is the Prandtl number, M = RBr, I — is the 


o 




Darcy number and Re is the Reynolds number. 

The corresponding non-dimensional boundary conditions 
are 

W = 0, 0 = 1 at r = 1 


56> 

= 0, — = 0 at r = 0 


dr 
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III. SOLUTION 


It is fairly unanimous that, the pumping action of the heart 
results in a pulsatile blood flow so that we ca represent the 
pressure gradient (pressure in the left ventricle) as 






and flow variables expresses as 
0(y,t) = 6 o (r)e' wt (9) 

w(y,t) = w () (r)e w, ( 10 ) 


Substituting Eqs. (8) - (10) into Eqs. (5) and (6) and solving 
the resultant equations subject to the boundary conditions in 
(7), we obtain 
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Here PI 2 — / co Pr Re— V 2 , 

p 2 = [M 2 + id) Re^(l + ) and 7 q (t) is the 


modified Bessel function of first kind of order zero. 

Thus the temperature distribution and the axial 
velocity are given by 
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The skin friction on the boundary of the tube is is 


given by 
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The rate of heat transfer coefficient in terms of 
Nusselt number on the boundary of the tube is 


Nu = 


d0 


dr 


_ Q/j (Q) 

= 4 («j 


e iG * (16) 


IV. DISCUSSION OF THE RESULTS 

Fig. 2 shows the effects of material parameter 
on W for M = 1, N = 1, p = 1 

co = 10 , A = 0.5, Pr = 0.7, Gr = 1, Re = 1 and 

t = 0.1. It is observed that, the axial velocity W increases 
at the axis of tube with increasing material parameter zl, , 

while it decreases near the tube wall with increasing . 

Effects of Hartmann number M on W for 

4=0.3 , N = 1, p = 1 , co = \0 , /l = 0.5, 

Pr = 0.7, Gr = 1, Re = 1 and t = 0.1 is shown in 

Fig. 3. It is found that the axial velocity W decreases with 
increasing Hartmann number M . 

Fig. 4 depicts the effects of Prandtl number Pr on 

W for A “ 0.3 

N = 1, p = 1 , oo = 10 , M= 1, 

zl = 0.5, Gr = 1, Re = 1 and ? = 0.1. It is noted that, the 
axial velocity W increases on increasing Prandtl 
number Pr . 
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Effects of radiation parameter JV on W for 

^=0.3 , p = 1 , co = 10 , Gr = 1, X = 0.5, 

M = 1, Pr = 0.7, Re = 1 and t = 0.1 is presented in 
Fig. 5. It is observed that, the axial velocity W increases with 
increasing N. 

Fig. 6 depicts the effects of Grashof number Gr on 

w for \ - 0.3 , TV = 1, p = 1, co —10 , M - 1, 
A = 0.5, Pr = 0.7, Re = 1 and t — 0.1. It is noted that, 

the axial velocity W increases at the axis of tube on 
increasing Grashof number Gr , while it decreases near the 
tube wall with increasing Gr . 

The effect of Reynolds number Re on W for 

TV = 1, \ =0.3, p = l,G) = 10, 1 = 0.5, M = 1, 

Pr = 0.7, Re = 1 and t — 0.1 is shown in Fig. 7. It is 
found that, the axial velocity W decreases with an increase in 
Reynolds number Re. 

Fig. 8 shows the effect of 1 on W for A [ = 0.3, 

TV = 1, M = 1, Pr = 0.7, Gr = 1, p = l,G) = 10, 

Re = l and t = o.i . It is observed that, the axial 

velocity W increases on increasing 1. 

Effect of oscillating parameter CO on W for 

M = 1 , p 0 = 1 , ^ = 0.5 , 1 = 0.5, Pr = 0.7, 
Gr = 1, TV = 1, Re = 1 and t = 0.1 is depicted in Fig. 

9. It is found that, the axial velocity W oscillates with CO . 

Fig. 10 shows the effect of Prandtl number Pr on 
6 for CO — 10, TV = 1, Re = 1 and t = 0.1. It is found 

that, the temperature 0 decreases with increasing Prandtl 
number Pr . 

Effect of radiation parameter TV on 0 for 
CO — 10, Pr — 0.7, Re = 1 and / = 0.1 is shown in 

Fig. 11. It is noted that, the temperature 0 increases with 
increasing N. 

Fig. 12 depicts the effect of Reynolds number Re 
on 0 for CO = 10, Pr = 0.7, TV = 1 and t — 0.1. It is 

observed that, the temperature 0 decreases with increasing 
Reynolds number Re. 

Effects of oscillating parameter CO on 0 
for Pr = 0.7, TV = 1, Re = 1 and / = 0.1 is depicted in 

Fig. 13. It is found that, the temperature 0 increases with 
increasing CO . 

Table-1 shows the effects of various physical 
parameters like , Gr , Re , X,M , Pr and TV on the 

skin friction T at the tube wall r — 1. It is found that, the 
skin friction increases for increasing the parameters , 

Gr , Re and A , while it decreases M , Pr and 
TV. r = 1. It is observed that, the Nusselt number increases 
with increasing Pr and Re , while it decreases with 
increasing N. 



Fig. 2. Effects of material parameter on XV for M = 1, 

p 0 = 1, co = 10, X = 0.5, Pr = 0.7, Gr = 1, 

N = 1, Re = 1 and t =0.1. 



Fig. 3. Effects of Hartmann number M on W for 

\ =0.5, p 0 =1, co — 10, A = 0.5, 

Pr = 0.7, Gr = l, TV = 1, Re = 1 and ? = 0.1. 
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Fig. 4. Effects of Prandtl number 1 on W for M = 1 , 

p 0 = 1, co = 10, A = 0.5, \ = 0.5, Gr = 1, 

N = l, Re = 1 and t =0.1. 



Fig. 5. Effects of radiation parameter on W for 

M =1, p 0 =1, co = 10,/l = 0.5, 

Pr = 0.7, Gr = 1, \ — 0.5, Re = 1 and t = 0.1. 



Fig. 6. Effects of Grashof number Gr on W for M = 1 , 

p 0 =1, co = 10, A = 0.5, Pr = 0.7, /l, = 0.5, 

N = 1, Re = 1 and t =0.1. 



Fig. 7. Effects of Reynolds number Re on W for M = 1, 

p 0 = 1, co —10, A = 0.5, Pr = 0.7, Gr = l, 

/V = 1, ^ = 0.5 and ? = 0.1. 



Fig. 8. Effects of pressure constant on W for M = 1 , 

p 0 = 1, co — 10,/^ = 0.5, Pr = 0.7, Gr = 1, 

TV = 1, Re = 1 and t = 0.1. 



Fig. 9. Effects of oscillating parameter CO on W for 

M =1, /? 0 =1, ^ =0.5,1 = 0.5, 

Pr = 0.7, Gr = 1, TV = 1, Re = 1 and ? = 0.1. 



Fig. 10. Effects of Prandtl number Pr on 0 for CO — 10, 

TV = 1, Re = 1 and t = 0.1. 
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r 

Fig. 11. Effects of radiation parameter N on 0 for 

co = \0,Pr = 0.7, Re = 1 and t = 0.1. 



Fig. 12. Effects of Reynolds number Re on 0 for 

CO = 10, Pr = 0.7, N = 1 and t = 0.1. 



Fig. 13. Effects of oscillating parameter CO on 0 

for Pv = 0.7, N = 1, Re = 1 and t = 0.1. 

Table - 1: The Skin friction T at the tube wall r = 1 



M 

Pr 

N 

Gr 

Re 

2 

T 

0.5 

1 

0.7 

1 

1 

1 

0.5 

-0.1118 

0.7 

1 

0.7 

1 

1 

1 

0.5 

-0.1072 

0.5 

2 

0.7 

1 

1 

1 

0.5 

-0.1189 

0.5 

1 

1 

1 

1 

1 

0.5 

-0.1308 

0.5 

1 

0.7 

2 

1 

1 

0.5 

-0.1609 

0.5 

1 

0.7 

1 

2 

1 

0.5 

-0.0809 

0.5 

1 

0.7 

1 

1 

2 

0.5 

-0.0814 

0.5 

1 

0.7 

1 

1 

1 

0.7 

-0.1072 


Table - 2: The Nusselt number Nu at the tube wall r = 1 


Pr 

N 

Re 

Nu 

0.7 

1 

1 

1.2113 

1 

1 

1 

1.6620 

0.7 

2 

1 

1.1200 

0.7 

1 

2 

2.0789 


V. CONCLUSIONS 

In view of these, we studied the effects of heat 
transfer and MHD on oscillatory flow of Jeffrey fluid in a 
circular tube with radiation. The expressions for the velocity 
field and temperature field are obtained analytically. It is 
found that the velocity W increases with increasing 

, Pr, TV, Gr, 69 and X , while it decreases with 

increasing M and Re; the temperature increases with 
increasing Pr ,N and CO, while it decreases with increasing 

Re ; the skin friction increases for increasing the parameters 
, Gv , Re and A , while it decreases M , Pr and N 

and the Nusselt number increases with increasing 
Pr and Re , while it decreases with increasing N . 
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